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The initial energy transfer in photosynthesis occurs between the light-harvesting pigments and 
on ultrafast timescales. We analyze the carotenoid to bacteriochlorophyll energy transfer in LH2 
Marichromatium purpuratum as well as in an artificial light-harvesting dyad system by using tran¬ 
sient grating and two-dimensional electronic spectroscopy with 10 fs time resolution. We find that 
Forster-type models reproduce the experimentally observed 60 fs transfer times, but overestimate 
coupling constants, which leads to a disagreement with both linear absorption and electronic 2D- 
spectra. We show that a vibronic model, which treats carotenoid vibrations on both electronic 
ground and excited state as part of the system’s Hamiltonian, reproduces all measured quantities. 
Importantly, the vibronic model presented here can explain the fast energy transfer rates with only 
moderate coupling constants, which are in agreement with structure based calculations. Counterin¬ 
tuitively, the vibrational levels on the carotenoid electronic ground state play a central role in the 
excited state population transfer to bacteriochlorophyll as the resonance between the donor-acceptor 
energy gap and vibrational ground state energies is the physical basis of the ultrafast energy transfer 
rates in these systems. 


I. INTRODUCTION 

Modern researchers have been fascinated by how effi¬ 
ciently solar photons are converted into chemical energy 
in photosynthesis. Initially, photons are absorbed in so- 
called light-harvesting complexes (LHCs). The following 
steps of the cascading energy transfer within LHCs are 
remarkably efficient: at low light intensities, 9 out of 10 
absorbed photons create a charge separated state in the 
reaction center as the basis for further energy conversion. 
m A thorough understanding of this process and its 
technological utilization will be a major contribution to 
a sustainable energy concept. Accordingly, this hope has 
motivated numerous studies, with the aim of finding and 
characterizing a bio-inspired, artificial light harvesting 
systems. [3j] Despite the efforts, artificial systems remain 
still less efficient and less stable than their natural coun¬ 
terparts. It is generally accpeted that a detailed picture 
of the involved electronic energy levels as well as the en¬ 
ergy dissipation pathways in both artificial and natural 
light harvesters is crucial in the realization of an efficient 
and durable bio-mimetic photosystem. 

While there is great structural diversity and flexi¬ 
bility in photosynthetic LHCs, [4j essentially only two 
molecular species serve as pigments, namely carotenoids 
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and (bacterio)chlorophylls, (B)Chls. Carotenoids and 
(B)Chls are fundamentally different in structure and 
complementing in physiological function. Carotenoids 
are linear molecules with varying endgroups. Their opti¬ 
cal properties are defined by the 7r-conjugated electronic 
states, extended along the polyene backbone as depicted 
in Fig. |T] The earliest suggested energy flow models used 
a three-level system with ground state So and excited 
states S 2 and Si [5|. The linear absorption spectra of 
carotenoids stem from transitions between the electronic 
ground state Sq and the “bright” electronic excited state 
S 2 , whereas the lowest-lying excited state .Si is “dark”, 
with the S 0 —> Si transition being one-photon forbid¬ 
den. The properties of Si are observed either by two- 
plioton absorption from So or by using nonlinear spec¬ 
troscopy, e.g. transient absorption (TA), where the ESA 

51 —1 S n is characteristically strong. @ The popula¬ 
tion flow rate from S 2 and Si strongly depends on the 
carotenoid’s chain length, with typical transfer times of 
less than 200 fs for solvated carotenoids (in vitro ).Besides 
quenching of harmful long-lived triplet states in (B)Chls 
and structural support of LHCs, the efficient transfer of 
the excitation energy from carotenoid to the lower lying 
(B)Chl states is the main functional role of carotenoids. 
The most efficient transfer route occurs from carotenoid’s 

5 2 state to the energetically closest (B)Chl-state, namely 
Q x . Transfer from Si to Q y is found in several an¬ 
tenna complexes and occurs on a longer timescale than 
S 2 —> Qx-M In vivo, i.e. as part of a LHC, the lifetime of 
S 2 decreases dramatically. [7J For example, in LHCII, the 
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main antenna system of higher plants and algae, the life¬ 
time of the S 2 state of the involved carotenoids is ~120 fs 
in vitro and only ^26 fs in vivo as measured by fluores¬ 
cence upconversion, which makes it one of the shortest 
in naturally occurring systems. Q The reason for this 
dramatic speed-up in carotenoid lifetime is the interac¬ 
tion with the closely neighboring (B)Chls. Chemically, 
the chromophore of (B)Chls is a nitrogen containing 
tetrapyrrole-macrocycle with one or two reduced pyrrole 
rings. According to the four orbital picture of porphyrins, 
0 the absorption spectrum of (B)Chls is described by at 
least two transitions in the so-called Soret band ( B x and 
By for HOMO - LUMO+1 and HOMO-1 - LUMO+1, re¬ 
spectively) and two Q-band transitions (HOMO - LUMO 
for Q y and HOMO-1 - LUMO for Q x ). Depending on 
molecular structure, the H-band is found in the blue to 
near-UV spectral region. The Q-bands stem from transi¬ 
tions with orthogonal transition dipole moments and are 
energetically degenerate for planar and symmetric por¬ 
phyrins such as metal-centered phthalocyanines. mm 
(B)Chls show a reduced symmetry in comparison, caus¬ 
ing a split of the Q-band. For BChl a in organic solvents, 
the Q x -band peaks near 600 nm, while the maximum 
of Q y is found just below 800 nm. Despite their struc¬ 
tural differences, carotenoids and (B)Chls fulfill comple¬ 
mentary tasks in photosynthetic antenna complexes. In 
vivo, excitonic coupling and the polar protein environ¬ 
ment shift the Q y -band of (B)Chls to the red, thus cov¬ 
ering the red edge of the visible solar spectrum. 

In this contribution we investigate an artificial dyad 
that mimics Rhodopseudomonas acidophila (LH2 Rps. 
ac.), and compare it to a naturally abundant LH2 of 
Marichromatium purpuratum (LH2 M. pur., formerly 
known as Chromatium purpuratum ). The main differ¬ 
ence between the artificial dyad and LH2 M. pur. is the 
smaller energy difference between S 2 and Q x in the latter. 
Based on this comparison, along with results reported in 
the literature for LH2 Rps. ac., we develop a novel energy 
transfer model for S 2 —> Q x energy transfer. We show 
that a transfer process mediated by vibronic resonances 
between S 2 and Q x presents a mechanistic alternative to 
conventional schemes incorporating only electronic cou¬ 
pling. 


II. MATERIALS AND METHODS 
A. Sample preparation 

The dyad-sample was provided by the group around 
Ana Moore and synthesized according to previously pub¬ 
lished procedures. 13 HI The dyad and its constituents 
were dissolved in spectroscopic grade toluene. 

Cells of Marichromatium purpuratum strain BN5500 
(also designated DSM1591 or 984) were grown anaero¬ 
bically in the light, harvested by centrifugation and sus¬ 
pended in 20 mM Tris-HCl pH 8.0. Then, upon the 
addition of DNase and MgCly, the cells were broken 


by passage through a French press. The photosynthetic 
membranes were collected from the broken-cell mixture 
by centrifugation and resuspended in the Tris-HCl so¬ 
lution to an optical density (OD) of 0.6 at 830 nm. 
The membranes were then solubilized by the addition 
of l%(v/v) of the detergent N, N-dimethyldodecylamine- 
N-oxide (LDAO). These solubilized membranes were sub¬ 
jected to sucrose density centrifugation with a step gradi¬ 
ent (as described previously by Brotosudarmo et al., 0) 
in order to fractionate the sample into LH2 and RC-LH1 
complexes. The LH2 fraction was then further purified by 
ion-exchange chromatography using diethylaminoethyl 
cellulose (DE52, Whatman). The bound protein was 
washed with 20 mM Tris-HCl pH 9.0 containing 0.15% 
dodecyl maltoside (DM) to exchange the LDAO deter¬ 
gent and then eluted by adding increasing concentrations 
of NaCl also containing 0.15% DM. Additional purifica¬ 
tion was achieved by gel-filtration chromatography us¬ 
ing a Superdex S-200 column (XK 16/100, GE Health¬ 
care) which had been equilibrated overnight with 20 mM 
Tris-HCl pH 9.0 plus 0.15% DM. 


B. Time resolved spectroscopy 

The two spectroscopic techniques employed in this ar¬ 
ticle, transient grating (TG) and two dimensional elec¬ 
tronic spectroscopy (2D-ES) are both four wave mixing 
techniques. In four wave mixing, three excitation pulses 
interact with the sample to create a non-linear polariza¬ 
tion of third order in the field, which acts as a source 
term for an emerging signal. In the experiment em¬ 
ployed here, all three excitation fields emerge along indi¬ 
vidual wave vectors, resulting in a non-linear signal direc¬ 
tion which is different from any of the incoming beams. 
This allows for a virtually background free signal. Details 
of the experimental layout are given elsewhere. fl5L ITgIJ 
Briefly, excitation pulses tunable throughout the visible 
spectral range are provided by a home-built non collinear 
optical parametric amplifier (NOPA), 0 pumped by a 
regenerative titanium-sapphire amplifier system (RegA 
9050, Coherent Inc.) at 200 kHz repetition rate. Pulse 
spectra were chosen to overlap with the investigated sam¬ 
ple’s absorption spectrum (see Fig.[T| and compressed 
down to a temporal FWHM of sub-11 fs in each case, de¬ 
termined by intensity autocorrelation. The pulses were 
attenuated by a neutral density filter to yield 8.5 nJ per 
excitation pulse at the sample. This corresponds to a 
fluence of less than 3.0 x 10 14 photons/cm 2 per pulse 
or 0.5% of excited molecules in the focal volume. The 
four wave mixing experiment used for both TG and 2D- 
ES relies on a passively phase stabilized setup with a 
transmission grating [18|, [l9| and has a temporal reso¬ 
lution of 5.3 as for t\, i.e. the delay between the first 
two pulses and 0.67fs for t 2 , the delay between second 
and third excitation pulse. A detailed description can 
be found in Milota et ai.[20] For TG measurements, t\, 
was kept at 0 fs while t 2 was scanned. The emerging 
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signal was spectrally resolved in W 3 by a grating-based 
spectrograph and recorded with a CCD camera. At a 
given delay time, spectra were not recorded on a shot- 
to-shot basis, but averaged over approximately 10 5 shots 
per spectrum. Sample handling was accomplished by a 
wire-guided drop jet [ 2 ll] with a flow rate of 20 ml/min 
and a film thickness of approx. 200 pm. The fact that 
the jet operates without the use of cell windows allows us 
to interpret 2D spectra and TG-signals even during pulse 
overlap (ti = f 2 = 0). All measurements were performed 
under ambient temperatures (295 K). 


damping constant A i E { 1 , 2 } and reorganization 
energy A BO,ki,i> * £ {1,2}. From the respective spectral 
density components | 22 j 
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line shape functions gBO,kl,i(t) °f th e bath can be calcu¬ 
lated by using the general formula for an arbitrary spec¬ 
tral density J(ui) , which reads 


C. Modeling 

1. Modeling of 2D-spectra 
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To facilitate the interpretation of the measured 20- 
spectra, results of model calculations were included in 
the analysis. It turned out that simulations including 
only carotenoid allowed most features of the measured 2D 
spectra to be explained. Therefore the model was based 
on the properties of the carotenoid component, whereas 
transfer from carotenoid to BChl only entered in terms 
of a contribution to lifetime broadening. With regards 
to electronic transitions it was assumed that the spec¬ 
tral profile of the pulses used in the experiment allows 
for resonant transitions between Sq and S 2 , and from Si 
to S n subsequent to population transfer from S 2 to Si 
. To account for vibrational effects in excitation pro¬ 
cesses involving the electronic states k and l, two modes 
were included as underdamped oscillators with reorgani¬ 
zation energies A uo,i> i £ { 1 , 2 } and vibrational frequen¬ 
cies ujuo »> i £ { 1 , 2 } in terms of line shape functions 


gUO,kl,i(t ) 
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The reorganization energies are connected to the Huang- 
Rhys factors Xuo,kl,i via A uo,kl,i = Xuo,kl,i^uo,i- 
Damping of the vibrations, which is known to appear in 
Si, was included by calculating the line shape function of 
the respective underdamped oscillator via double integra¬ 
tion of a correlation function } 22 |, where the real part of 
the correlation function was approximated by including 
terms up to a selected order n of a series characterized 
by the Matsubara frequencies v n = 27Tk ^ Tn . 

Furthermore, during the evolution in electronic states 
k and l the influence of the bath was described by two 
Brownian oscillator spectral density components with 


The sum of all line shape function components was in¬ 
cluded in the response functions for the calculation of 
2D-spectra. The response functions contain the transi¬ 
tion frequencies uju between the electronic states k and 
l and rate constants of transfer processes. Besides the 
rate ks 2 ^Si f° r transfer from S 2 to Si also a transfer 
rate ks 2 ^Q n: from S 2 to an electronic state Q x of the 
Chlorophyll component of the LH2 complex plays a role 
in lifetime broadening effects. 

The rate of lifetime broadening in S 2 consists of the 
sum i(fc 5 2 _ > 5 1 +fcs 2 ->-Qx) an d is denoted as Tg 2 . For life¬ 
time broadening in Si a rate was introduced. The 
response functions can be distinguished by their rephas¬ 
ing or nonrephasing properties and assigned to excita¬ 
tion processes of stimulated emission (SE), ground state 
bleaching (GSB) or excited state absorption (ESA) type, 
where in the case of SE and GSB only the electronic tran¬ 
sition dipole moment gs 0 S 2 enters, while for ESA com¬ 
ponents with population transfer the transition dipole 
moment gs 1 S Tl a ls° appears. 

The response functions were formulated in analogy to 
[,23| and (24j and are given in the Appendix [I] The sum of 
all rephasing and nonrephasing response function compo- 
nents yields the third-order response Sjj ( 73 , 72 ,ti) and 

/o\ 

^nr( t - 3 > r 2 , Tl ), respectively, which reflect the excitation 
dynamics in the limit of infinitely short pulses. The in¬ 
fluence of the finite pulse width pulses was taken into 
account in analogy to Ref. [25;j, i.e. by convolution of 
the calculated signal with the pulses. By integration 
over the time intervals n, T 2 and 73 between the elec¬ 
tronic transitions, the resulting third-order polarization 
is transformed into a dependence on the pulse delays t \, 
t 2 and £ 3 . The rephasing part of the 2D-spectrum is ob¬ 
tained via 
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the nonrephasing part results from 


02D,NR(wi,t2,W3) 




dt 3 exp(iwiti) 


x exp(iu 3 t 3 )P$ ) R (t 1 ,t 2 ,t 3 ). (5) 


2. Calculation of Energy Transfer Rate in a Vibronic Model 

To simulate the energy transfer rate from the 
carotenoid S 2 state to the Q x of the BChl in LH2 and the 
dyad we employ a hetero-dimer model. This approach is 
natural for the dyad showing now signs of aggregation be¬ 
tween the cliromophores. For LH2 the monomeric treat¬ 
ment of Q x and S 2 is also justified, given that there are 
no signs of excitonic splitting for these bands. Moreover 
the carotenoid predominantly interacts with one of its 
neighboring B850 BChls with the next nearest neighbor 
coupling showing only half the interaction strength. j26j 

To model the carotenoid molecule we use an effec¬ 
tive two state model which involves two electronic states 
(So and S 2 ) and one effective vibrational mode which 
replaces the two fast vibrational modes with frequen¬ 
cies fli ft 1150 cm” 1 and Q 2 ~ 1500 cm -1 known to 
be present the carotenoid and it simulates the progres¬ 
sion of states |0)i5oo|0)ii5o, |0) 15001 1) 1150 ? |l)i50o|0)n5o, 
J 0 )i 50 o| 2 }ii 50 ,11) 1500 1 1 ) 1150 , |2) 1500 |0)n 50 ) ■ • • which orig¬ 
inates from the two modes. We can group these levels ac¬ 
cording to vibrational energies equal approximately av¬ 
erage of two mode values to 0, ft 1325 cm -1 , ft 2650 
cm” 1 , etc. (counted from the zero point energy). Be¬ 
cause the spectral lines corresponding to the transitions 
to the states listed above are broadened, they act ef¬ 
fectively as one mode with a roughly average frequency. 
This carotenoid mode (referred to as primary mode fur¬ 
ther on in this paper) is coupled to a bath of harmonic os¬ 
cillators |27H29j (a secondary modes or secondary bath) 
which provides damping and dephasing. The coupling 
between the effective mode and the carotenoid electronic 
transition provides the optical dephasing leading to the 
broadening of the absorption spectra. In order to cal¬ 
culate correctly the energy transfer rate, we parametrize 
the fast mode and its coupling to the bath of harmonic 
oscillators in such a way that it reproduces absorption 
spectra. This effective model of the bath enables us to 
limit the bath description to one harmonic mode whose 
dynamics is treated explicitly by a master equation ac¬ 
cording to m ■ Th this way the mode effectively provides 
a correct bath spectral density for the energy transfer 
between the carotenoid and BChl while allowing for a 
non-perturbative treatment. Fit of the absorption spec¬ 
tra of the carotenoid in the dyad yields the carotenoid 
effective vibrational frequency w e ff = 1390 cm” 1 , HR 
factor s e ff = 1.3, the reorganization energy of the sec¬ 
ondary bath X b = 670 cm” 1 and the bath correlation time 


Tb = 30 fs. For the LH2 M. pur. we obtained w e g = 1453 
cm” 1 , s e ff = 1.21, A b = 1242 cm” 1 and r b = 54 fs. 

The Q x states of the BChl and purpurin are treated by 
the same model. The representative mode on the BChl 
and purpurin molecules are chosen so that they fit ab¬ 
sorption spectra including the significant vibrational side 
bands of the molecules. The side band for the purpurin 
can be easily located in the monomeric absorption spec¬ 
trum given in Fig. [T] By fitting the purpurin monomeric 
absorption spectrum we get the effective mode frequency 
Upurp = 1246 cm” 1 , HR factor s purp = 0.61, reorgani¬ 
zation energy X purp = 600 cm” 1 and bath correlation 
time of T purp = 47 fs. BChl absorption spectrum ex¬ 
hibits a much weaker vibrational side band and leads to 
a smaller HR factor of sbcm = 0.22. Other parameters 
are as follows: toschi = 1284, Xbcm = 1911 cm” 1 and 
tBC hl = 100 fs. The effective vibrational modes repre¬ 
sent the overdamped part of the spectral density found 
in BChls (see e.g. Ref. [30]) by the term discussed in 
Section HV B 21 and the side band by its first replica. The 
strength of this replica is given by the HR factor. 

The parameters obtained by the fit of the monomers 
are slightly readjusted in a fit of the aggregate spectra, ei¬ 
ther the purpurine-carotenoid dyad or the LH2 complex. 
Energy transfer rates are then estimated with no further 
adjustment of the parameters. To estimate the energy 
transfer rate, the population dynamics of the system is 
calculated by the master equation from Ref. [29| with 
the initial population corresponding to the laser pulse 
excitation of the carotenoid. The total time dependent 
population of the Q x state which is initially zero is calcu¬ 
lated and the population transfer time is defined as the 
time ts 2 ^q x when the population of Q x state is equal to 
1/e . 


III. RESULTS 
A. Absorption 

Fig. Ha) shows the absorption spectrum of the dyad 
(black) and its constituents, i.e. the carotenoid donor 
(green), and the purpurin acceptor (red). 

Details on the dyad’s structure and dynamics are given 
in detail elsewhere. dill mi It is a donor-acceptor sys¬ 
tem, consisting of a /3-carotene derivative (donor, green in 
the molecular structure in FigOJ a ))> and a tetrapyrrole- 
macrocycle as an acceptor, referred to as purpurin (red). 
The dyad-carotenoid and purpurin are linked by an 
amide bond, providing structural rigidity but only par¬ 
tial conjugation between donor and acceptor. The dyad 
was synthesized to mimic LH2 Rps. ac., whose struc¬ 
ture is depicted in Fig.QJb), as taken from the Papiz et 
o/,'.32| LH2 Rps. ac. consists of nine pairs of pigment- 
protein subunits (a //3-subunits), forming two concentric 
cylinders, to which carotenoids and BChl a molecules are 
bound non-covalently. As can be seen in Fig. Hb), the 
BChls are arranged in two ring-like structures (red and 
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Figure 1. (a) Absorption spectra and molecular structure of 
the caroteno-purpurin dyad (black line) and its constituting 
carotenoid (donor, green) and porphyrinic acceptor in red. (b) 
Structural model and absorption spectrum (c) of LH2 Rps. 
ac. (d) Absorption spectra at different temperatures of LH2 
M. pur., along with the molecular structure of the bound 
carotenoid (okenone). 


yellow) with different numbers of BChls (18 for red, 9 
for yellow) and ring diameters, leading to a weakly and 
strongly interacting set of BChls (yellow and red, respec¬ 
tively). As a result, the excitonic part of the absorp¬ 
tion spectrum of LH2 Rps. ac., shown in Fig.QJc), ex¬ 
hibits two peaks, one near the monomeric Q^-transition 
at 800 nm and a red-shifted peak near 850 nm. The red 
and yellow sets of BChls in Fig. [TJb) are therefore re¬ 
ferred to as B850 and B800, respectively. The carotenoid 
is in van der Waals contact (<3.5'i^ce) with both the 
B800 and the B850 ring. [32l [iisl] The dyad lacks such ex¬ 
citonic bands, as its structure is monomeric and there are 
no interacting BChl-moieties. An interesting difference 
between the absorption spectrum of the dyad and LH2 
Rps. ac. arises in the carotenoid region of the spectrum. 
The vibronic structure of the monomeric carotenoid’s ab¬ 
sorption spectrum, shown in green in Fig. (Ha), is more 
pronounced than in the dyad, shown in black. This can 
be understood by invoking vibronic coupling between 
donor and acceptor, as will be addressed in detail in sec¬ 
tion IIVI 

The absorption spectrum of LH2 M. pur. (Fig.QJd)) 
exhibits substantially less vibronic modulation in 
carotenoid region of the spectrum. This is attributed to 
the polar endgroup of the bound carotenoid (okenone), 
interacting strongly with the polar in vivo protein 
environment. As the bright SVstate carries ionic 
character,[34j its energy relative to So is strongly affected 
by the interaction with the protein and shifts to the red in 
comparison the carotenoid (rhodopin glucoside) of LH2 
Rps. ac. The BChl a molecules in LH2 M. pur. don’t 
show similar bathochromic shifts, which leads to a greatly 
enhanced overlap of the carotenoid- and Q^-band in this 
system. The excitonic region of the spectrum of LH2 M. 


pur. exhibits a less pronounced red shift with respect to 
the monomeric transition as compared to LH2 Rps. ac. 
Recent preliminary X-ray structure analysis [35jj explains 
this by a reduced number of a//3-subunits (eight for LH2 
M. pur. and nine in the case of LH2 Rps. ac.) 


B. Transient Grating 

Fig. UKa) shows the emission frequency (W 3 ) resolved 
transient grating (TG) signal of LH2 M. pur. The sig¬ 
nal shows a broad frequency range, spanning the entire 
spectrum of the excitation pulse, shown as an orange line 
in Fig. [21)b), in comparison to the absorption spectrum 
of LH2 M. pur. depicted as a grey area. The double- 
peaked structure of the excitation spectrum explains the 
overall shape of the TG-signal. The dashed lines mark 
frequencies characteristic for S 2 (dashed green line) and 
Q x (dashed dark red line). Figures [U(c) and (d) show 
the respective transients. For both detection frequen¬ 
cies, we observe a fast decaying peak around t 2 = 0 fs, 
followed by a slower mono-exponential decay with a de¬ 
cay constant of 330 ± 50 fs according to least square fit 
analysis. The signal does not decay to zero within the 
1000 fs time window of the experiment. The observed 
decay behavior is, within error margins, identical for 
both detection frequencies. TG probes the sum of the 
absorptive and dispersive part of the induced third or¬ 
der signal, while pump-probe selects only the absorp¬ 
tive part, meaning that the retrieved time constants can¬ 
not be directly compared to pump-probe measurements. 
[371 . -37] Vibrational dynamics however will yield simi¬ 
lar frequencies between the two techniques (see. i.e. vi¬ 
brational dynamics in /3-carotene probed by pump-probe 
[38j j and by TG ®). The vibrational response mani¬ 
fests as an oscillatory signal, superimposed on a slowly 
varying background. After subtraction of the latter and 
Fourier transformation of the remaining signal along t 2 
for every detection frequency, the (u >2 , w 3 ) frequency 
map in Fig. [51(e) is retrieved. A cut along a ^-specific 
detection frequency (green dashed line) yields the well- 
known carotenoid frequencies as depicted in Fig.^f ) near 
1000 cm -1 (Methyl-rocking motion), 1160 cm -1 (carbon 
single bond stretching) and 1530 cm -1 (carbon double 
bond stretching), and are in agreement with Resonance 
Raman spectra measured for okenone. |40j The same 
set of frequencies is observed when detecting at the Q x - 
band, as shown in Fig.[5Kg). This is interesting because 
BChl-specific modes near 1600 cm -1 or 700 cm -1 as know 
from resonance Raman measurements EH are not found. 
The absence of BChl vibrational signatures is readily ex¬ 
plained by vibrational-electronic coupling in BChl, which 
is small in comparison to carotenoids. [42j 
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Figure 2. (a) Detection frequency ( 013 ) resolved transient 

grating signal of LH2 M. pur. (b) Absorption spectrum in 
gray in comparison to the spectrum of the excitation pulse 
in orange. Temporal cuts at ^-specific detection frequencies 
specific for S 2 and Q x in (c) and (d), respectively, (e) Fourier- 
transform spectra along 12 , dispersed along LJ 3 . (f) and (g) 
show cuts through (e) for S 2 - and (jU-specific detection fre¬ 
quencies, respectively. 


C. 2D electronic spectroscopy 

As mentioned in section eh 2D-ES and TG are related 
in the sense that they are both four wave mixing tech¬ 
niques. The difference lies in the treatment of the first 
inter-pulse delay t\, which is kept at zero in TG, but 
scanned and Fourier-transformed (ti — > oj\) in 2D-ES. 
The emerging (wi, W 3 )- 2 D plots, recorded at a fixed sec¬ 
ond inter-pulse delay t 2 , allow for a correlation of excita¬ 
tion (wi) and emission frequencies ( 013 ). Fig- 0 compares 
electronic 2D spectra at £2 = 0 fs for the dyad and LH2 
M. pur. 

At t 2 = 0 fs, electronic population is given no time to 
relax to lower lying states, within experimental time reso¬ 
lution. Both spectra in Fig. [3] show strong positive (red) 
peaks along the diagonal ( 011 = 013 ). Comparison to the 
absorption spectra, given in the lower panels in Fig. EKa) 
and (b), show that the strong diagonal peaks stem 
from the carotenoid’s So —>■ S 2 transition. The dyad’s 
spectrum in Fig.[5Ja) shows two well resolved peaks 
corresponding to the Q^-transition at approximately 
17200 cm ” 1 and a carotenoid peak near 19000 cm” 1 . The 
off-diagonal or cross peaks (wi 7 ^ 3 ) with negative ampli¬ 
tude in Fig. [5]stem from ESA-patliways. For carotenoids, 
ESA from S 2 in the visible has been observed. j43j' , [44] 
The ESA of the dyad’s porphyrinic acceptor was de¬ 
scribed previously as broad and featureless. While the 
combination of donor and acceptor ESA explains the neg¬ 
ative features in the dyad’s 2D-spectrum, we note that 
at t 2 = 0 fs, negative signals do not necessarily have to 
be related to ESA, as they already occur in a two level 
system at this waiting time. J45i] Positive off-diagonal 
peaks can stem from vibrational energy levels. |16l l46H48j 
stimulated emission (SE) or coupling between two elec¬ 
tronic oscillators. [451] The latter would be the most inter- 


(a) Dyad, 0 fs 



1- 1 -1-•-r 

17 18 19 

co. / 1000 cm' 1 

1 


(b) LH2 M. pur. , 0 fs 



Figure 3. 2D electronic spectra of the dyad (a) and LH2 M. 
pur., each at f 2 = 0 fs. Positive (negative) signals are drawn 
a red (blue) full lines at 5% steps. Each figure is normalized 
to its respective maximum. The dashed grey lines mark the 
nodal line at zero intensity. The lower panels in each graph 
show the overlap between absorption (dark red line) and ex¬ 
citation spectrum (green line). 


esting scenario, as electronic coupling between donor and 
acceptor would suggest an excitonic energy level struc¬ 
ture and the associated transfer mechanisms.[49] Upon 
visual inspection, the location of the positive cross peaks 
in the dyad’s 2D signal suggests electronic coupling as 
a likely explanation. The stronger of the two positive 
cross peaks, found below the diagonal (wi > 0 J 3 ), peaks 
at a value of W 3 that coincides with the absorption fre¬ 
quency of the Q^-band. There is also a corresponding 
peak above the diagonal (wi < 013 ), which peaks blue de¬ 
tuned with respect to Q x and with weaker intensity than 
its counterpart below the diagonal. Such differences in 
intensity may however be attributed to the ultrafast pop¬ 
ulation transfer from S 2 to Q x , the onset of which takes 
place within the experimental time resolution of Ilfs. 
The blue detuned maximum of the upper cross peak can 
be attributed to finite pulse width or effects of chirp in 
the excitation pulses. (50) Inspection of the 2D signal of 
LH2 M. pur disproves the assumption of electronic cou¬ 
pling. The major difference between the dyad and LH2 
M. pur is the position of the ^-transition. While the 
dyad’s carotenoid exhibits maximum absorption above 
19000 cm” 1 , the carotenoid in LH2 M. pur, okenone, has 
a carbonyl-containing endgroup as shown in Fig.njd), 
shifting the 0-0 transition of okenone in vivo down to 
17570 cm” 1 , while the energetic position of the Q x - 
band remains unchanged between dyad and LH2 M. pur. 
Hence, the energy difference between S 2 and Q x should 
be greatly reduced in LH2 M. pur, and the correspond¬ 
ing cross peaks should be shifted, leading to a large and 
fairly featureless electronic 2D-signal. In contrast to this 
expectation, the lower cross peak in LH2 M. pur shows 
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roughly the same energy difference from the diagonal of 
approximately 1500 cm” 1 . This strongly suggest a vibra¬ 
tional origin, given that the prominent C=C stretching 
mode of carotenoids is found at this energy. 

The peak above the diagonal in LH2 M. pur is harder 
to explain. The ESA-signal in Fig.|5Jb) is much weaker 
than for the dyad case, which is readily explained by a 
different ESA-transition energy from £2 for okenone in 
vivo. (36] We can thus observe a weak diagonal peak at 
wi = W 3 ~16480cm, with an excitation frequency 
corresponding to the upper cross peak. This energetic 
position excludes coupling to Q x as a possible explana¬ 
tion, peaking at 17035 cm ” 1 in LH2 M. pur. The physical 
origin of the upper cross peak in Fig.|5Jb) becomes more 
obvious when examining the peak shape evolution along 
t 2 , as shown in Fig.|4] 

The upper diagonal peak is discernible up to 30 fs, but 
drops to zero by 80 fs. This behavior is unexpected for 
an electronic coupling peak, whose amplitude should only 
be a function of coupling constant J, but not of waiting 
time t 2 - The rapid decay of the upper cross peak rather 
suggests that it is an artifact, related to the shape of the 
excitation pulses. The excitation spectra are clearly non- 
Gaussian (see lower panel of Fig.|3](b)), making limited 
pulse width or chirp related effects [50j the most likely 
explanation for the upper diagonal cross peak for LH2 
M. pur. 

Another interesting feature is observed at (aq = - 
17500cm” 1 , ui 3 = 16480cm” 1 ), i.e. an energetic position 
indicative of S 2 —> Q x energy transfer pathways. The 
peak builds up rapidly within the first 30 fs, which cor¬ 
responds well with the 55 fs transfer time reported for 
this process. (SSI After 50 fs, negative ESA-signals domi¬ 
nate the spectra. This is explained by the intramolecular 
S 2 —>■ £1 energy transfer, occurring at a 95 fs build up 
rate as reported previously. (36j We support this assign¬ 
ment by calculating the electronic 2D-spectra only for 
okenone, the carotenoid of LH2 M. pur. Details of the 
calculations can be found in section HTC. 


IV. DISCUSSION 

A. Carotenoid-Chlorophyll Interaction 

Spectroscopic experiments reveal the transfer from Car 
£2 and £1 states to Chls and BChls as fast and relatively 
efficient despite the rather short £2 state life time of the 
monomeric carotenoid (at least in comparison with Chls 
and BChls). Both Q x and £2 states are optically al¬ 
lowed, and they can therefore interact by the resonance 
coupling mechanism. Because of the proximity of the two 
molecules in the LH2 structure, one cannot exclude an 
electron exchange (Dexter) mechanism. However, quan¬ 
tum chemical calculations suggest Coulomb interaction 
to be dominant here (52(. In this case, theoretical de¬ 
scription the carotenoid-BChl interaction and modeling 
of the corresponding excitation energy dynamics remain 


within the well understood framework of excitonic de¬ 
scription of photosynthetic antennae fdljl [53 ]. 

In the standard framework, we assume the electronic 
coupling element between the two electronic states to 
be independent of the vibrational DOF of the molecu¬ 
lar system. Electronic excitation is transferred between 
the two molecules due to an interplay of the resonance 
coupling and the fast fluctuations of the molecules’ re¬ 
spective energy gaps. The energy gap fluctuations are 
caused by their interaction with the nuclear degrees of 
freedom (DOF) which involve both the intra- as well as 
intermolecular vibrations. These are considered a ther¬ 
modynamic heat bath characterized by certain temper¬ 
ature and spectral density. Two regimes of the relative 
coupling strengths are usually distinguished. The delo¬ 
calized regime in which the system-bath interaction is 
weak in comparison with the resonance interaction be¬ 
tween the molecules, and the localized regime in which 
the system-bath coupling is considered strong. The terms 
localized and delocalized refer to the effective electronic 
eigenstates of the interacting system. For the delocal¬ 
ized regime the coupling between the two molecules is 
strong enough to sustain correlation between the elec¬ 
tronic excitations on different molecules so that delocal¬ 
ized excitons are formed. In the localized regime (also 
referred to as Forster regime here), the system-bath in¬ 
teractions induce energy gap fluctuations which prevent 
any such correlation, and the excitations appear to be 
localized on the molecules. In fact, even in relatively 
weakly excitonically coupled systems one can often de¬ 
tect signs of delocalization (see e.g. Ref. [54[), because 
the actual apparent eigenstates are always somewhere in 
between the two theoretical limits. For photosynthetic 
aggregates formed by BChl and Chi molecules, where the 
protein bath can be to some extent assumed as unstruc¬ 
tured, modern simulation methods such as the Hierar¬ 
chical Equations of Motion (HEOM) enable us to deter¬ 
mine the excited state dynamics beyond the two limits 
mentioned above. jUj However, even in the case of ag¬ 
gregates composed of BChl or Chi molecules only, some 
clearly visible spectroscopic features can result from the 
involvement of pronounced vibrational modes. j56i f50| 

In the discussion that follows, we want to argue that 
the involvement of the vibrational DOF of the carotenoid 
molecules is crucial for the energy transfer dynamics be¬ 
tween the £2 and the Q x states. We will compare two 
theoretical approaches to the calculation of the energy 
transfer rates, both motivated by the observed features 
of the spectra and known spectroscopic behavior of the 
system. We will argue that both support the crucial role 
of the fast vibrational modes, in particular the carotenoid 
ground state vibration levels. 

The Hamiltonian of the relevant part of the system 
(state Qx of the BChl and the £2 state of the carotenoid) 
can be written in the following way 

H = H,< q b J + Hg) + (e Qx + AVq.(</«■>)) \Q x ){Qx\ 



Figure 4. Electronic 2D spectra of LH2 M. pur., measured at indicated values of £ 2 - Line coloring follows the same conventions 
as described in Fig. [3]. 


+ (es 2 + AVs 2 (g (S2) )) |S 2 )(S 2 | 

+ J r (|S 2 )(Q x | + |Q x )(5 2 |). (6) 

Here, cq x and es 2 are the optical electronic transi¬ 
tion energies of the BClil Q x state and the carotenoid 
S '2 state, respectively (including the reorganization en¬ 
ergy of the bath), H^ and are the Hamil¬ 

tonian operators of the nuclear DOF on the BChl 
and carotenoid, respectively. These two sets on 
nuclear DOF described by macroscopic sets of co¬ 
ordinates gWD = {qi Chl , ■ ■ ■, qn Chl } an d q^ S2 ' > = 

{qi ar ,..., q^ ar }, where n is a macroscopically large 
number, form the bath interacting through operators 
A VQ x (q^ x ' 1 ) and AVs 2 (q^ S2 ' > ) with the electronic tran¬ 
sitions on the BChl and carotenoid, respectively. The 
electronic coupling element J describes the resonance 
interaction between the collective excited state | Q x ) = 
\ e Q x )\gCar) in which the BChl is excited to the state 
Q x and the carotenoid is in its electronic ground state 


and the collective excited state |5 2 ) = \gBChi)\e-s 2 ) in 
which the BChl is in its electronic ground state and the 
carotenoid is excited to its excited state 5 2 . To com¬ 
plete the Hamiltonian we should also to include the state 
I/) = |eQ x )|e S2 ) in which both Car and BChl are excited. 
However it will be shown later that the features related 
to this state are negligible as the resonance coupling will 
be shown to be weak. In the experimental 2D spectra, 
shown in Fig. U one can easily identify the fast rising 
ESA component, which is usually assigned to the strong 
absorption from the 5i states of the carotenoid. Con¬ 
tribution of this state is taken into account by including 
a separate partial Hamiltonian and S 2 —> S 1 relaxation 
rate into the response functions in Appendix |T] All other 
states with no spectroscopic contributions are taken into 
account in the simulations by introducing decay rates for 
the bright states. For the subsequent discussion it will 
be useful to introduce a bare electronic (single exciton) 
Hamiltonian which reads as 

H e i = £q x \Qx)(Qx\ + es 2 |‘S 2 )(S' 2 | 
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+ J(\S 2 )(Q x \ + \Q x )(S 2 \). (7) 

This Hamiltonian represents the basic three level (ground 
state energy was chosen to be zero) structure of the prob¬ 
lem. 


B. Energy Transfer Rates 

As discussed above, the interplay of energy gap fluc¬ 
tuations and the resonance coupling energy determines 
which theoretical limit applies. The usual discussion of 
the relative strengths of the two interactions is made 
based on the relative values of the reorganization energy 
A which describes the magnitude of the energy gap fluc¬ 
tuations and the value of the resonance coupling J [60']. 
Our system is characterized by reorganization energies 
A q x of the BChl and As 2 of the Car which should be 
compared to the resonance coupling energy J. For the 
validity of the weak resonance coupling limit it is usually 
required that A > \J\. In our case, assuming J « 100 
cm -1 (see Ref. |26l l6lj l confirms the validity of this 
regime because the total A s 2 is in the order of thousands 
of cm -1 as it involves two fast vibrational modes with 
Huang-Rhys factors close to 1 and vibrational frequen¬ 
cies above 1000 cm -1 as shown in Fig. [2] The Q y state 
of BChl in the protein environment of the LH2 complex 
exhibits A Q y ~ 100 cm -1 (see e.g. Ref. [62j l and we 
can assume a similar value A q x ~ 100 cm' 1 as for the 
Q x state. Here, the reorganization energy of the chro- 
mophore is only comparable, not larger, than the res¬ 
onance coupling element. We note that conventionally 
the role of temperature is completely neglected in this 
discussion. The magnitude of the energy gap correlation 
function, which provides a true measure of the energy 
gap fluctuations, depends on temperature and the depen¬ 
dency is linear in the high temperature limit (in the case 
of Brownian oscillator model, see Ref. It would 

therefore be more suitable to compare the value of A ksT 
and | J| 2 . At room temperature, ksT = 207 cm' 1 and 
A Q x ksT > |J| 2 . Yet another problem arises when consid¬ 
ering the role of the relative energy gap Ae = es 2 — cq x 
of the Q x and S 2 states. In an ideal dimer, the delo¬ 
calization is determined by the so-called mixing angle 
9 = b arctan p^j.j64l| The delocalization decreases with 
the increasing Ae. The interplay of the relative energy 
gap, temperature, reorganization energy and resonance 
energy is a subject of numerous theoretical studies in the 
context of photosynthesis. However, no simple formula 
taking into account the influences of all these parameters 
on the degree of delocalization exists. 

After the discussion of the parameters of the 
carotenoid-BChl dimer, we have many reasons to be¬ 
lieve that the interaction of the S 2 and Q x states can be 
treated in the weak coupling limit. The issues is however 
more complicated because of the nature of the system- 
bath interaction in the carotenoid. Most of its reorgani¬ 
zation energy can be assigned to the fast intramolecular 


vibrational modes. The spacing between the vibrational 
levels and the relative electronic energy gap As are such 
that the Q x state can be viewed as effectively interact¬ 
ing with the nearest vibrational levels and not with the 
carotenoid electronic state as a whole. This is a situ¬ 
ation similar to the one studied in Refs. [56j, j5f| H5 
where special resonances between vibrational and elec¬ 
tronic energy gaps lead to pronounced effects despite the 
small Huang-Rhys factors of the BChls. Clearly, treating 
the carotenoid vibrations as a thermodynamic bath is a 
questionable approximation. 

In the following subsections we will briefly discuss the 
carotenoid-BChl energy transfer rates based on the weak 
coupling Forster theory, and compare it with an explicit 
treatment of the carotenoid vibrations in the Hamilto¬ 
nian, along with Refs. [56| i59f. The energy level struc¬ 
ture characteristic for the two situations is presented in 
Fig. El 


1. Weak Resonance Coupling Limit 


In the weak resonance coupling limit, the parame¬ 
ter J of the Hamiltonian, Eq. ©, is assumed small 
and perturbation theory to the second order can be 
performed. For two transitions (donor and acceptor), 
the transition rate is given by the Fermi golden rule 
as Ka<-d = ^\J\ 2 5(tu — £a)- It is possible to view the 
donor molecule as a set of many available de-excitation 
transition (the molecule is excited electronically and vi- 
brationally) while the acceptor molecule can be viewed 
as a set of transitions ready to be excited. In a realis¬ 
tic case it is therefore necessary to integrate the Fermi 
formula over all acceptor and donor transitions 

Ka^d = ~r\J \ 2 
n 


OO OO 

J de D J de A f ahs (e A )h{e D )S(e D -e A ), 


( 8 ) 


where / a b s and /a correspond to the normalized distri¬ 
bution of the transition energies available for excitation 
on the acceptor and the normalized distribution of the 
transition energies available for the deexcitation on the 
donor. The Forster rate can be derived directly from 
Eq. © by noticing that the acceptor absorption spec¬ 
trum is related to the distribution of transition frequen¬ 
cies /abs(w) = fabs{kio) as a(w) ~ w/ a bs(w) and similarly 
for the donor fluorescence spectrum cr(w) ~ w 3 /a(w) (see 
e.g. [53l l&jjl or directly by perturbation theory with re¬ 
spect to J and a cumulant expansion for the bath DOF 
as in Ref. [49]. In both cases there is a direct relation 
between the rate and the overlap of absorption spectrum 
of the acceptor with the emission spectrum of the donor 

K -o J J l 2 / duaA(u)a D (uj)u >- 4 

Ai ~ D h 2 J (lujaA(w)u }- 1 f cIuj<jd(oj)uj - 3 
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Figure 5. Interaction schemes of the electronic and vibra¬ 
tional levels in (a) localized (Forster) and (b) delocalized (vi- 
bronic) regime. In the case (a) the transition energy from the 
non-equilibrated state of the carotenoid excited state vibra¬ 
tions has to be resonant with the Q x transition energy. The 
relative energy gap Ae between the S 2 and Q x transitions 
requires involvement of more than one vibrational quanta 
of the carotenoid. The vibronic case (b) also involves exci- 
tonic mixing of the vibrational levels of the carotenoid with 
the Q x transition. The collective single excited states of the 
carotenoid-purpurin or carotenoid-BChl dimer involve states 
in which the Q x transition is excited and simultaneously the 
electronic ground state vibrations of the carotenoid are ex¬ 
cited. The states can excitonically mix with the carotenoid 
electronically excited state with a different number of vibra¬ 
tional quanta. 


It is important to note that the Forster theory does not 
in general depend on the dipole-dipole coupling approx¬ 
imation although the specific dependence of the rate in 
this approximation is often used. From the point of view 
of the present manuscript, the coupling energy J can be 
calculated by any type of advanced methods of quantum 
chemistry which takes into account all the details of the 


electronic structure of the excited states of the interact¬ 
ing molecules [65]. From the theoretical point of view, 
the Forster rate can be derived from Hamiltonian, Eq. 
©, under various assumptions. Usually one assumes 
relaxed state of the environment and the intramolecu¬ 
lar vibrational modes. However, Forster theory is very 
versatile and allows many generalizations [65l 1661 ], most 
importantly it is still valid for the case that the excited 
state of the donor is unrelaxed ^3j. Eq. (J5J is then still 
valid, but it does not result in Eq. ©• Forster theory 
therefore also allows for a direct experimental estimation 
of the resonance coupling if the fluorescence spectrum of 
the donor and the absorption spectrum of the acceptor 
are known, even in case that the fluorescence does not 
start from a relaxed state. For the present dyad, the 
coupling J was determined to be J = 240 cm -1 , within 
the framework of Eq. © . [T3 | 

Fig. Eta) describes a typical situation for the short 
time dynamics of the S 2 ~Q X energy transfer. The higher 
vibrationally excited levels of the S 2 are populated, and 
in order to achieve resonance with the Q x transition, even 
higher vibrational levels on the carotenoid ground state 
are excited after deexcitation of the donor (S 2 ). The en¬ 
ergy transfer is enabled by the fact that the energy cor¬ 
responding to the Ae is accepted by the ground state vi¬ 
brational levels of the carotenoid molecule as depicted in 
Fig. ]5](a). The same statement can be rephrased in terms 
of the spectral overlap. Because the rate depends on the 
overlap of the donor emission and the acceptor absorp¬ 
tion spectra, the effective broadening of the carotenoid 
spectrum due to the fast vibrational modes increases sig¬ 
nificantly the energy transfer rate. Forster theory thus 
assigns the electronic ground state vibrational levels of 
the carotenoid a crucial importance in the energy trans¬ 
fer process. 


2. Vibronic Coupling 

The main conclusion of the preceding section is that 
the vibrational modes of the carotenoid are the main 
driving element of the energy transfer from S 2 to Q x . 
Forster theory shows that these modes provide resonance 
for the energy transfer and that this resonance involves 
multiple vibrational quanta. As one can learn from Fig. 
(Ha), absorption lineshape undergoes a slight change in 
the carotenoid region suggesting weak delocalization be¬ 
tween Q x and S 2 states, at least for dyad. It is there¬ 
fore natural to include the important vibrational modes 
explicitly into the Hamiltonian, and to treat only the 
remaining overdamped modes as a bath. This leads di- 
rectjy to the model of vibronic excitons (see Refs. [561 - 
|58jj|67|)- Such a model automatically treats the selected 
vibrational modes non-perturbatively. 

We split the nuclear part of the carotenoid Hamilto¬ 
nian in Eq. © into some selected modes denoted by 
coordinates q n and referred to later in the text as pri¬ 
mary modes. The rest of the modes which we denote by 
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q'(S 2 ) ( see g ec |i y A I for the notation) will be referred to 
as secondary modes or secondary bath. We have 

H S B 2 = H sl + tfvib + ffvib-b • (10) 

where is the Hamiltonian of the secondary bath and 
the Hamiltonian of the primary modes and their interac¬ 
tion with the secondary bath read as 

= ail 

n 

and 

H v ib-b = ~ ^ ( Qn ^ ( Cfc Qk- (12) 

n k 

Here, u> n are the frequencies of the primary modes to 

be treated explicitly, and qk are the coordinates from the 
set qd s ^) coupled to the primary modes by some coupling 
constants In the excited state of the carotenoid, the 
bath Hamiltonian is now split in the following way 

+ A V s ,tf 8a) ) = H$> + A Vs 2 (q' {S2) ) 

+ E “W (Pl + (Qn - dn) 2 ) - E ZkQk- 

n n k 

(13) 

Here, we assume that the same secondary bath which 
drives the ground state vibrations into the equilibrium, 
drives also the excited state vibrations to their corre¬ 
sponding equilibrium. Correspondingly, the energy gap 
operator reads as 

A V S2 (q {S2) ) = A Vs 2 (q' {S2) ) 

- E tuVndnqn + E dn E ^ qk ■ ( 14 ) 

n n k 

We assumed that the spectral densities describing the in¬ 
teraction of the selected modes with the rest of the bath 
are the same. The energy gap operator, Eq. (HD, is 
composed of three parts, A Vs 2 (qd S2 ' > ) representing the 
overdamped bath, the part representing the direct influ¬ 
ence of the selected modes on the energy gap fluctuations 

Hel —vib = ^ ( fvdjndnqm (15) 

n 

and the part describing the influence of the overdamped 
bath on the electronic energy gap functions via the cou¬ 
pling with the selected modes. The strength of the latter 
part of the interaction Hamiltonian is given by the shift 
d n of the excited state vibrational potential with respect 
to the ground state potential. The essence of the vibronic 

model is to include H e i_ v ib to the Hamiltonian which we 

treat explicitly and only an effective energy gap operator 

AHjfV^)) = A Vs^q'^) + E d "E^ ( 10 ) 

n k 


is treated by perturbation theory. It leads to the bath 
correlation function 

C(t) = (AVs 2 (q' {S 2 ) ;t)AVs 2 (q' (S 2 ) -,0)) + ci el \t), (17) 

where the last term on the right hand side corresponds to 
the direct contribution of the secondary bath to the de¬ 
phasing. The damping of the nth primary modes is gov¬ 
erned by the bath correlation function originating from 
the term, Eq. (ED, 

Ci n \t) = Y / \^\ 2 (Qk(t)qk), (18) 

k 

while the direct contribution of the secondary bath to the 
dephasing is expressed by the same correlation function 
scaled by the dimensionless shift of the effective vibra¬ 
tional mode 

^b el) (*) = Ei d «i 2c b" ) (*)- (19) 

n 

The calculation of the energy transfer dynamics in the 
vibronic model thus starts with a diagonalization of the 
Hamiltonian 

H s = Hvib + H el _ vih \S 2 )(S 2 \ + H e i. (20) 

The diagonalization results in eigenstates which combine 
the electronic and vibrational states into the mixed vi¬ 
bronic states. In a second step the rates of population 
transfer and coherence dephasing are calculated for the 
set of eigenstates by second order perturbation theory 
and cumulant expansion [27] using correlation function, 
Eq. (fTSl) . for the damping of the selected mode and the 
correlation function on ly to simulate electronic 

dephasing. This means that we assume the correlation 
function C(t), Eq. (fT71) . to be composed of the contri¬ 
bution of the secondary bath only. This choice decreases 
the number of degrees of freedom in fitting and effectively 
ties strength of the coupling between the primary mode 
and the secondary bath to the strengths of the interac¬ 
tion between the electronic transition and the primary 
mode. Because only effective spectral density is impor¬ 
tant for the calculation of the transition rates, we believe 
that the advantage in limiting the number of parameters 
outweighs the restriction put on our model. In princi¬ 
ple one could estimate the overall rate of energy transfer 
directly by weighting the calculated rates between indi¬ 
vidual levels. We have however chosen the more practical 
approach of calculating the transfer time as described in 
Section m~C2l 

Because of the rather large vibrational frequency and 
moderate coupling, the Q x transition can be viewed as 
interacting only with the energetically nearest transition. 
This situation is depicted in Fig. [Ub). For the sake of 
clarity of the graphical presentation, we depict a situa¬ 
tion of perfect resonance between the relative energy gap 
Ae and the vibrational frequency and show energy level 
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splitting of the energetically nearest levels. Fig. [5}o de¬ 
picts the relevant carotenoid and BChl energy levels as 
two separate entities (in gray) and in terms of collective 
states of the carotenoid-BChl dimer in the case of J = 0. 
Conceptually, the ground state vibrational levels of the 
carotenoid contribute to an excited collective state if the 
Q x state is excited. This is a consequence of the fact 
in a singly excited aggregate state all molecules except 
one are in their ground state (carotenoid in this case). 
The ground state vibrational levels of the carotenoid 
thus interact resonantly with the vibrational levels of the 
carotenoid in the electronically excited state S 2 ■ This in¬ 
teraction corresponds to the energy transfer from S 2 to 
Q x with a simultaneous deposition of the excess energy to 
the ground state vibrations of the carotenoid. The com¬ 
plete picture of the energy levels of the system is much 
more complicated because the excitonic mixing involves 
all available states. The details of this mixing, which 
are included in the basis transformation from the local¬ 
ized states to the energetic eigenstates, lead to discernible 
changes in the lineshape of the carotenoid molecule. As 
will be demonstrated in the following sections, the vi- 
bronic model provides an explanation for the observed 
spectral changes as shown in Fig. Oja). Also, the de¬ 
tailed treatment of the interaction between carotenoid 
vibrational levels and the Q x state leads to the exper¬ 
imentally observed energy transfer rates with a signifi¬ 
cantly weaker resonance coupling than predicted by Eq. 
©• 


C. Simulation of Absorption Spectra and Transfer 
Rates in Vibronic Model 


The vibronic model discussed in the previous section 
was implemented as described in Section III C 21 The 
absorption spectra of monomeric carotenoid was fitted 
with a model Hamiltonian involving one effective high 
frequency mode to obtain starting values of the effective 
mode parameters for the subsequent fitting of the ab¬ 
sorption spectra of the dyad and the LH2 complex. We 
found the frequency of the vibrational mode in the dyad 
to be w e ff = 1390 cm” 1 and its Huang-Rhys factor was 
determined to be s e ff = 1-3. All monomeric parameters 
are discussed in Section III C 21 

In Fig. [S]we summarize our fitting efforts for the dyad. 
We chose the dyad because the same solvent can be used 
for the dyad and its components avoiding solvent related 
bathochromic shifts. The absorption spectra of LH2 M. 
pur. and LH2 Rps. ac. were equally well reproduced 
with the employed vibronic model. Fig. H a ) presents 
the result of summation of the dyad’s monomer spec¬ 
tra. This is the spectrum corresponding to the weak 
resonance coupling case (Forster regime) and it is pre¬ 
sented to highlight the changes between monomeric and 
dyad spectra. In Fig. H>, we present a result of the fit¬ 
ting in which the Huang-Rhys factor of the carotenoid, 
the Q x transition frequency, the relative energy gap Ae 


between S 2 and Q x and the resonance coupling J were 
allowed to change. The fitting achieves good agreement 
with the experimental dyad spectrum, most importantly 
in the carotenoid related part. There the relative hight 
of the absorption maxima of different vibrational peaks 
changes between the monomer and the dyad. The es¬ 
timated coupling is J = —119 cm -1 while s e g = 1.29, 
Ae = 2400 cm -1 and 6 q x = 17150 cm -1 . The minor 
change with respect to the parameters of the monomers, 
such as the change of the transition energy to the excited 
state, have to be assigned to the chemical change upon 
formation of the dvad.(l2| 

The same fitting procedure with the fixed resonance 
coupling J estimated from measured fluorescence spec¬ 
trum by Forster theory, J = 240 cm -1 [l2|, results in a 
less satisfactory agreement as can be seen in Fig. He). 
In order to demonstrate that lineshape changes in all 
cases are not only due to the change of the monomeric 
Huang-Rhys factor of the carotenoid molecule, we show 
the spectrum of the dyad with the parameters resulting 
from the same parameters but with J = 0 cm -1 in Fig. 
Etd). Interestingly, the best fit is achieved with a change 
of the Huang-Rhys factor of less than 1 % and the change 
in the line shape can be therefore mostly assigned to the 
resonance coupling. 
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Figure 6. Absorption spectra of the dyad. Panel (a) shows 
the sum of the monomeric spectra, which corresponds to the 
weak coupling case (no changes to the lineshape except of 
lifetime broadening); (b) presents the best fit of the dyad 
spectrum with the Q x transition energy, relative energy gap 
Ae, carotenoid HR factors and the resonance coupling free to 
be adjusted. Panel (c) shows the best fit with the resonance 
coupling fixed at J = 240 cm -1 and panel (d) presents the 
sum of the monomers with the parameters corresponding to 
the best fit from panel (b), but with J = 0 cm -1 . 


The simulation of the absorption spectrum has resulted 
in a value of J which is by 50 % smaller than the one esti¬ 
mated from the experiment, suggesting that higher value 
is not compatible with the spectral changes with respect 
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to the monomers seen in the dyad. The obvious question 
is whether this lower value for the resonance coupling is 
compatible with the measured ultrafast rates. To this end 
we perform calculations of the S 2 to Q x transfer dynam¬ 
ics to estimate the transfer time introduced in Section 
III C 21 for various values of J and Ae and keeping other 
parameters of the model fixed to those from the best fit 
from Fig. [^b). The plot of transfer time as a function 
of J and Ae is presented in Fig. [7] The calculations 
yield a transfer time of ts 2 ^q x = 55 fs for the values of 
the best fit which is in good agreement with the exper¬ 
imental value of 40 fs.|12j] Moreover the 2D plot in Fig. 
[7] and the corresponding cut at the value of J = —119 
cm _1 reveals a periodic modulation of the transfer time 
with the period corresponding to the effective spacing 
between vibrational levels. The transfer time decreases 
monotonically with increasing value of J, as shown in 
the lower panes of Fig. [7] representing cuts along the ex¬ 
perimental vales of Ae. The periodic modulation of the 
transfer time is in agreement with the resonant involve¬ 
ment of the carotenoid vibrational modes in the energy 
transfer. As in the Forster case this implies the involve¬ 
ment of the carotenoid vibrational levels from the elec¬ 
tronic ground state. Conservation of energy during the 
energy transfer requires the excess energy to be deposited 
to the vibrational energy of the carotenoid. 


D. Simulation of the 2D Spectra 

In order to assign the observed features of the ex¬ 
periment 2D-spectrum we performed simulations of the 
carotenoid 2D spectra. This is a first step towards a 
complete simulation of the interacting Car-BChl system, 
which allows us to assign Car only features and possibly 
identify 2D spectral components originating in Q x state. 
Such a full calculation will be a subject of our future 
work. 

Calculating the 2D spectrum in a spectral region cov¬ 
ering the edge of the studied molecule is an extremely 
difficult task given the fact that a good fit on the outer 
region of (even the absorption) spectrum depends on 
many rarely studied details of the chromophore inter¬ 
action with its environment. For instance, assumptions 
about the Gaussian disorder of the electronic transitions 
energies and typical estimates for bath spectral densi¬ 
ties are well suited for explaning the features around the 
maximum of absorption and do not cover well the situ¬ 
ation at the edges of the spectra. On the other hand, 
introducing more fitting freedom by assuming arbitrarily 
some other types of spectral densities and disorder dis¬ 
tributions seems not to yield more theoretical insight, as 
the new parameters cannot be fixed by a limited set of 
experiments. Below we therefore describe a calculation of 
the Car 2D-spectra in which the Car molecule should well 
represent the Car component of the LH2. Correspond¬ 
ingly a limited fitting to the experiment 2D-spectra from 
Fig. Q] is done, keeping in mind that the two systems 
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Figure 7. Transfer time as a function of J and Ae for the dyad 
(a) and for LH2 M. pur. (b). At the experimental value of 
Ae and the values of resonance coupling determined by fitting 
absorption spectrum, transfer time for the dyad and LH2 M. 
pur. are 55 fs and 40 fs, respectively, in agreement with the 
values reported in pump probe measurements. 


which are compared in such a fitting are different. Large 
part of the Car parameters are therefore motivated by 
the values from literature. 


1. Carotenoid only Spectra 

From the linear absorption spectrum of the LH2 com¬ 
plex the electronic excitation energy ws 2 s 0 = 19770 cm -1 
was obtained (see Fig. [[]). Because of the line shape func¬ 
tion approach, which relies on second-order cumulant ex¬ 
pansion with the electronic ground state as the reference 
state, this value corresponds to the vertical transition 
energy, i.e. the sum of the electronic excitation energy 
and the reorganization energies of the spectral density 
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components. The transfer rates were chosen under the 
assumptions 1 /A;g-,= 95 fs and 1 /ks 2 ^Q x = 55 fs in 
agreement with [68|. Because of the slow transfer from 
Si to other electronic states, the lifetime broadening con¬ 
stant Tg-, was taken as zero. Even though some spectral 
density parameters for okenone are available from the 
literature (68j, the respective parameters for /3-Carotene 
have been reported in more detail. They can at least 
provide some orientation for the choice of the parame¬ 
ters to model the okenone component of the investigated 
LH2 complex. According to [63 the frequencies of the 
included vibrational modes are ujjjo.i = 1150 cm -1 and 
u>uo,2 = 1520 cm -1 . The Huang-Rhys factors of the vi¬ 
brational modes in S2 were assumed as Xuo,S 2 S 2 ,i = 0.25 
and Xuo,S 2 S 2 ,2 = 0.5. While these values are somewhat 
smaller than the ones for /3-Carotene given in Ref. [ 69j , 
their ratio is similar. According to the tendency reported 
in the literature, the Huang-Rhys factors in Si were as¬ 
sumed to be larger than those in S 2 by a scaling factor 
chosen as 1.5. To describe the damping of the vibra¬ 
tional modes in Si with a time constant of 50 fs, ten 
Matsubara terms were taken into account in the calcu¬ 
lation of the lineshape function. For one of the Brow¬ 
nian oscillator spectral density components a fast decay 
with damping constant 1 /Asoi,S 2 S2 = 30 fs and a reor¬ 
ganization energy Xbo,S 2 S 2 ,i = 300 cm -1 in agreement 
with Ref. @ was assumed. To reproduce the inhomo¬ 
geneous broadening effects in the absorption spectrum of 
the LH2 complex, the second Brownian oscillator spec¬ 
tral density component with very slow decay constant 
of 20 ps (73 was included in the model. For a reorga¬ 
nization energy of \bo,S 2 S 2 ,2 = 3000 cm -1 comparable 
inhomogeneous broadening effects as in a measured lin¬ 
ear absorption spectrum could be obtained and features 
in the measured 2D-spectra could be qualitatively repro¬ 
duced. In Ref. 0 it has been reported that coherence 
gets lost during this population transfer process, so that 
the fluctuations in S 2 and Si can be considered as uncor¬ 
related. Accordingly, gs 2 s 1 and gs 2 s n were taken as zero. 
The Huang-Rhys factors for excitation from Si to S n are 
known to be much smaller than the ones of the excitation 
from So to S 2 . In Ref.j68| a value of 0.1 is given, which 
was assumed for both modes in our calculation. 

According to our interpretation, this value of the 
Huang-Rhys factor is given with respect to So, as the 
vibrational energy for the transition between Si and S n 
reported in Ref. [68| cannot be explained by such a small 
Huang-Rhys factor. However, the much larger Huang- 
Rhys factors of Si with respect to Sq allows for a larger 
vibrational energy than expected from the value of 0.1 
for the Huang-Rhys factor. 

Regarding the Brownian oscillator modes, it has been 
reported in 0 that for double excitation from Si the 
reorganization energy corresponds to a fraction of only 
0.5 of the corresponding reorganization energy in Si from 
the electronic ground state, which was taken as the same 
as for S 2 in our description. 

The vertical transition energy between Si and S„ was 


chosen as 0 Js 1 s n = 15400 cm 1 in agreement with Ref. 
6 §|. In the prefactors of the ESA response functions 
2 

= 1.5 entered, following the values given in 

■ other carotenoid derivatives. The finite pulse 
widths were taken into account in terms of Gaussian 
profiles determined from a fit of the measured local os¬ 
cillator spectrum . In this way the central frequency 
ujq = 16275 cm -1 and a FWHM of 1612 cm -1 were ob¬ 
tained. As the fitted curve corresponds to the squared 
pulse profile, multiplication of the FWHM with a factor 
of \J2 was required to obtain the FWHM of the single 
pulses with a value of 2280 cm -1 . 

In the real part of the sum of all calculated response 
function contributions after convolution with the pulses 
an intensive, positive-valued diagonal peak in the re¬ 
gion of (w 1 = —17500 cm - 1 ,W 3 = 17500 cm -1 ) appears 
at all population times up to 100 fs, starting from 
t 2 = Ofs. This peak mainly stems from the rephas¬ 
ing GSB contribution. Furthermore, below the diago¬ 
nal at (wi = —17500 cm - 1 ,W 3 = 16000 cm -1 )a crosspeak 
from the rephasing GSB contribution is found, which 
stems from vibrational effects. This initially oval peak 
becomes butterfly-shaped with a nodal line separating 
positive an negative region at £2 = 10 fs, whereas at 
t 2 = 20 fs it becomes completely positive again. At 
t 2 = 30 fs a modified shape of the vibrational crosspeak 
below the diagonal from the rephasing GSB contribution 
in combination with a rise of of the rephasing SE con¬ 
tribution and the negative-valued rephasing ESA contri¬ 
bution leads to a localization of the crosspeak close to 
(uq = —18000 cm - 1 ,W 3 = 15500 cm -1 ). From f 2 = 50 fs 
the ESA contribution starts to obscure the latter. At 
t 2 = 80 fs the negative ESA peak covers a broad ener¬ 
getic range below the diagonal. This tendency becomes 
even more pronounced at £2 = 100 fs and t 2 = 200 fs, 
where the ESA peak pushes the diagonal peak almost 
completely above the diagonal. 

Different from the measured 2D-spectra, the upper di¬ 
agonal peak at t 2 = 0 to 30 fs is missing in the calcu¬ 
lated 2D-spectra shown in Fig. [5] most likely due to the 
fact that the spectral phase of the pulses was considered 
as flat in the calculations. When comparing the spec¬ 
tra from experiment (Fig.[]| and simulation (Fig. [5]) at 
t 2 = 30 fs, it becomes apparent that the experimentally 
obtained cross peak, which we identified as an indica¬ 
tion of S 2 —> Q x energy transfer, is missing in simula¬ 
tions. Considering that the simulations only incorporate 
the okenone-part of the spectrum, this missing cross peak 
confirms our assignment to energy transfer to Q x . 2D-ES 
elucidates why carotenoid to bacteriochlorophyll energy 
transfer pathways are difficult to analyze in pump-probe 
spectroscopy. 0 0 The according cross peak coin¬ 
cides with a recurring negative carotenoid feature (see 
t 2 = 30 fs and £2 = 10 fs spectrum in Fig.0. Additionally, 
S 2 —» Q x energy transfer is overlaid with ESA from Si, 
making it only observable within the ultrashort lifetime 
of S 2 , i.e. only during 55 fs in LH2 M. pur. 
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Figure 8. Simulated electronic 2D spectra of okenone, the carotenoid in LH2 M. pur. The employed values of t 2 are indicated 
in the panels. Line coloring follows the same conventions as described in Fig. [3] 


V. CONCLUSIONS AND OUTLOOK 

In this work we have shown that a vibronic cou¬ 
pling mechanism describes several central aspects of the 
carotenoid to BChl energy transfer dynamics, yielding 
moderate coupling constants J close to 100 cm” 1 * * . This 
value is in excellent agreement with structure-based cal¬ 
culations and reproduces experimental absorption spec¬ 
tra well. The low values of J obtained from the vi¬ 
bronic model suffice for explaining the experimentally 
observed ultrafast transfer rates. This is explained by 
vibronic resonances, which have the potential of dramat¬ 
ically speeding up energy transfer. [4Sj, |g3, [zij Hence, 
the vibronic coupling mechanism circumvents the short¬ 
comings of Forster theory and explains static as well as 
dynamical properties of natural light harvesters. 

The investigated carotenoid-BChl system is an extreme 
case of a heterodimer, meaning that donor (carotenoid) 
and acceptor (BChl) differ in transition energies and 
HR-factor. Future studies will test the vibronic cou¬ 
pling mechanism on energy transferring dimers where 


the spectroscopic properties are more alike, such as pery- 
lene bisimide dyads. (73 In such systems, both donor and 
acceptor show strong electron-phonon coupling, which 
should make vibronic effects in energy transfer dynamics 
even more pronounced. Another highly promising class 
of systems is bulk heterojunction solar cells, where vi¬ 
bronic coupling between donor and acceptor was recently 
suggested to be the mechanism behind ultrafast electron 
transfer. [?6| 
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I. RESPONSE FUNCTIONS 

In this Appendix, the response functions used to calculated the Car 2D-spectrum are specified. The rephasing GSB 
term with population of S 2 and a decay due intramolecular transfer to S± is given as 

R3g,s 2 ( T 3,T2,Ti) = |/hs 0 s 2 1 4 exp(iw,s 2 s 0 Ti — *ws 2 s 0 t 3 ) exp (— Tg 2 Ti -Ts 2 t 3 ) 

ex P(—5s 2 s 2 ( r i) + 9*s 2 s 2 ( t 2 ) - 9s 2 s 2 (t 3 ) 

~9s 2 S 2 ( t 1 + r 2) — 9s 2 S 2 ( r 2 + r 3) + 9s 2 S 2 ( t 1 + t 2 + T3)). (1) 


The rephasing SE component with decay of the S 2 component due to intramolecular transfer to Si reads 
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R2g,S 2 (T3,T 2 ,Ti) = \bs 0 S 2 | 4 exp(*Ws 2 s 0 n - iUs 2 S 0 T 3 ) exp (-Ts 2 Ti - r S2 T 3 ) 

ex P(-5s 2 s 2 ( T i) + 9s 2 s 2 {t 2 ) - g*s 2 s 2 (,T 3 ) 

-g*s 2 s 2 ( r i + T 2 ) - gs 2 s 2 (r 2 + r 3 ) + ffs 2 s 2 (n + r 2 + r 3 )) 

exp(-fcs 2 ^s 1 r 2 ). (2) 

For the rephasing ESA component with population transfer from S 2 to S\ and subsequent excitation to a higher 
excited state one obtains 

R if,s 2 ^s 1 ( T 3,T2,Ti) = - |/xs 0 s 2 | 1 2 3 \ns 1 s n \ 2 exp(?'ws 2 s 0 n - iws„s 1 r 3 )exp(-r S2 n - r Sl r 3 ) 
ex P(—5 s 2 s 2 ( t i) - 9s 2 sA T A - ffSiS^Ts) 

+9s 2 s i ( r i + r 2) + gs 2 s 1 (T 2 + t 3 ) — gs 2 s! ( r i + r 2 + r 3 )) 
exp (gs n s 2 (t 2 ) + 2 ff5nSl (r 3 ) - g* SnS2 (n + r 2 ) 

-5s n S 2 (t 2 + r 3 ) + ffs nS2 (n + r 2 + r 3 ) - g Sn s n (r 3 )) 

PT 2 

/ dr exp(—fcs^.-s.gj.r) 

Jo 

exp(2 i^(gs 2 s l (r 2 - r) - gSiSi (r 2 - r) + psiSi (r 2 - r + r 3 ) - ^s 2 Si (r 2 - r + r 3 ))) 
exp(2i3(gs rl s 1 (r 2 - r) - gs n S 2 (r 2 - r) + gs n s 2 {j 2 - t + r 3 ) - gs n sA T 2 - r + r 3 ))). 

The nonrephasing response functions of GSB-, SE- and ESA-types read 

R ig,s 2 (r 3 , r 2 , ri) = |^ So s 2 1 4 exp(-*ws 2 s 0 n - ius 2 s 0 r 3 ) exp (-r S2 n - r S2 r 3 ) 

exp( gs 2 s 2 (r i) gs 2 s 2 A 2 ) gs 2 s 2 (r 3 ) 

+gs 2 s 2 (n + r 2 ) + 3 S 2 S 2 (r 2 + r 3 ) - gs 2 s 2 (n + r 2 + r 3 )), (3) 

d?i s ,s 2 (r 3 ,r 2 ,ri) = |^ So s 2 1 4 exp(-*ws 2 s 0 n - ius 2 s 0 r 3 ) exp (-F^n - r S2 r 3 ) 

exp(-ss 2 s 2 (ri) - g*s 2 s 2 {T 2 ) ~ g*s 2 sA T s) 

+5s 2 s 2 (n + r 2 ) + gs 2 s 2 ( r 2 + r 3 ) - ffs 2 s 2 (n + r 2 + r 3 )) 

exp(-fcs 2 ^ Sl r 2 ) (4) 

and 

R 2 f,s 2 ^yS! ( T 3> r 2 ,n) = - |ms 0 s 2 | 2 lMs 1 s Il | 2 exp(-*ws 2 s 0 ri - iu} SnSl T 3 ) exp (-r S2 n - r Sl r 3 ) 
exp(-5s 2 s 2 (ri) +ffs 2Sl (r 2 ) - gs 1 S l {r 3 ) 

~9s 2 Si (n + r 2 ) - ffs 2iSi (r 2 + r 3 ) + 3 s 2Sl (n + r 2 + r 3 )) 
ex p(— 5 s n g 2 (r 2 ) + 2 ff5nSl (r 3 ) + g Sn s 2 (n + r 2 ) 

+5s n s 2 ( t 2 + r 3 ) - ffs n s 2 (n + r 2 + r 3 ) - (r 3 )) 

PT 2 

/ drexp(-fcs 2 ^ Sl r) 

Jo 

exp(2 i^(gs 2 S! (r 2 - r) - gSiSi (r 2 - t) + psiSi (r 2 - r + r 3 ) - ^s 2 5i (r 2 - r + r 3 ))) 
exp(2i3(gs n s 1 (r 2 - r) - gs n s 2 (r 2 - t) + gs„s 2 (r 2 - r + r 3 ) - gs n sA T 2 ~t + r 3 )))- 
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